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Abstract 

Gauged linear sigma models with C m -valued scalar fields and gauge group U(l) d , d < m, 
have soliton solutions of Bogomol'nyi type if a suitably chosen potential for the scalar fields 
is also included in the Lagrangian. Here such models are studied on (2 + 1) -dimensional 
Minkowski space. If the dynamics of the gauge fields is governed by a Maxwell term the 
appropriate potential is a sum of generalised Higgs potentials known as Fayet-Iliopoulos 
D-terms. Many interesting topological solitons of Bogomol'nyi type arise in models of 
this kind, including various types of vortices (e.g. Nielsen-Olesen, semilocal and supercon- 
ducting vortices) as well as, in certain limits, textures (e.g. CP m_1 textures and gauged 
CP m ~ textures). This is explained and general results about the spectrum of topologi- 
cal defects both for broken and partially broken gauge symmetry are proven. When the 
dynamics of the gauge fields is governed by a Chern- Simons term instead of a Maxwell 
term a different scalar potential is required for the theory to be of Bogomol'nyi type. The 
general form of that potential is given and a particular example is discussed. 
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1. Introduction 

In the study of topological solitons, field theories of Bogomol'nyi type occupy a spe- 
cial place. Mathematically such field theories are characterised by the fact that their 
soliton solutions may be obtained by solving first order differential equations (called the 
Bogomol'nyi equations) instead of the second order Euler-Lagrange equations. A further 
characteristic property is the existence of static multisoliton solutions made up of arbi- 
trarily placed single solitons. Physically one may think of this property as reflecting the 
absence of static forces between well-separated single solitons. It follows that, in theories of 
Bogomol'nyi type, the space of static multisoliton solutions is larger and more interesting 
than in generic field theories with topological solitons. Defining the moduli spaces to 
be the spaces of static soliton solutions of soliton number N, one finds for example that 
in theories of Bogomol'nyi type the dimension of Mn increases linearly with N. Starting 
with the proposal by Manton |1J in the context of magnetic monopoles, moduli spaces 
have become a powerful tool for investigating the physics of several interacting solitons in 
theories of Bogomol'nyi type . This program is particularly advanced in the case of vortices 
in the abelian Higgs model. Here moduli space techniques have been used successfully to 
study the interacting dynamics of several vortices || and, more recently, the statistical 
mechanics of vortices ||. 

The purpose of this paper is to explore a certain class of (2+1) -dimensional field the- 
ories with topological solitons of Bogomol'nyi type which provide a natural generalisation 
of the abelian Higgs model. In two spatial dimensions, it is useful to distinguish two types 
of topological solitons: vortices and textures (depending on the context the latter are also 
often called lumps or baby-Skyrmions). Generally we may distinguish the two by the origin 
of their topological stability. If the stability is due to a non-trivial first homotopy group 
of the vacuum manifold we speak of vortices, if it is due to a non-trivial second homotopy 
group of the vacuum manifold we speak of textures (though care has to be taken in in- 
terpreting these statements, as first pointed out in ref. |§ and as we shall also see later). 
More specifically, vortices arise in field theories with charged scalar fields when a U(l) 
gauge symmetry is spontaneously broken, whereas textures typically arise in non-linear 
sigma models. Typical and much studied examples of such textures are the lumps in CP™ 1 
sigma models. 

Remarkably, both vortices and textures can be studied in a unified manner in the 
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framework of gauged linear sigma models (GLSM's) with a judiciously chosen Higgs po- 
tential. The special potentials which are required here are known to supersymmetry the- 
orists as Fayet-Iliopoulos D-terms, and GLSM's with these potentials have been studied 
intensively in recent years in the context of topological sigma models and string theory, 
see for example refs. HH- It has also been emphasised in that context that, in suit- 
able limits, GLSM's with Fayet-Iliopoulos D-terms become non-linear sigma models whose 
target spaces are certain special manifolds called toric varieties (of which the projective 
spaces CP m are particular examples). Thus there exists a unifying framework for studying 
textures and vortices of Bogomol'nyi type, with a host of beautiful mathematical results, 
whose implications for the study of topological defects in (2+1) dimensional field theories 
do not appear to have been fully exploited. 

However, while there is significant overlap between the questions studied by string 
theorists and those which one might ask in the context of (2+1) -dimensional field theories, 
there are also important differences. While in the context of topological sigma models and 
string theory GLSM's have typically been studied on compact Riemann surfaces, the field 
theorist would naturally study the models on flat Minkowski space. Thus, in the (2+1)- 
dimensional case static fields are defined on (non-compact) R 2 , and it is natural also to 
consider time-dependent fields. This leads to additional questions. On a non-compact 
domain the convergence of integrals needs to be checked carefully. This results in extra 
constraints which affect the spectrum of Bogomol'nyi solitons. The inclusion of the time 
coordinate allows one to study time-dependent fields and to consider gauge fields whose 
dynamics is governed by a Chern-Simons term. The latter is physically interesting because 
topological solitons in such theories display anyonic statistics. 

This paper begins with two general sections in which GLSM's are introduced and gen- 
eral results about the spectrum of Bogomol'nyi solitons are derived. Most of the algebraic 
manipulations in these sections are standard, but the parts of the analysis which deal with 
the non-compactness of R 2 appear to be new. In sect. 4 we show how the non-linear CP 1 
sigma model can be understood in a very precise way as a limiting case of the semilocal 
vortex model in the context of GLSM's, and we explain how this limiting procedure can 
be generalised. Our contribution in this section is mainly an expository one, linking work 
done in the study of topological defects with the mathematical framework of toric varieties. 
In sect. 5 we study in detail a family of models with two complex scalar fields and gauge 
group U(l) x U(l). We show that the topological solitons in this family include, for various 
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parameter values and in various limits, superconducting vortices and gauged CP lumps. 
The latter are solitons in the gauged 0(3) sigma model introduced in [(?[. The present 
paper developed out of an attempt better to understand the mathematical structure and 
physical interpretation of that model, and we will see that the framework of GLSM's with 
Fayet-Iliopoulos D-terms provides a satisfactory understanding of both. In sect. 6 we write 
down a general Lagrangian for GLSM's of Bogomol'nyi type with Chern-Simons terms for 
the gauge fields. Again we illustrate the general results in a particular model with two 
complex scalar fields. Finally, the purpose of sect. 7 is to draw together the rather diverse 
viewpoints which enter this paper and to highlight some open questions. 

2. Gauged linear sigma models 

Perhaps the most natural way to introduce GLSM's with Fayet-Iliopoulos D terms 
without invoking supersymmetry is to recall the manipulations that lead to the estab- 
lishment of the Bogomol'nyi equations for the abelian Higgs model. The basic fields in 
the abelian Higgs model are a complex scalar and a U{1) gauge field, but as we shall see 
presently Bogomol'nyi equations can still be established, mutatis mutandis, if we consider 
instead a C m -valued scalar field w and gauge any subgroup of the maximal torus of the 
unitary group U(m) acting on w in the fundamental representation. The maximal torus of 
U(m) is m-dimensional and a choice of generators t a , a = 1, ...m, defines an isomorphism 
between it and U(l) m . Concretely the generators are diagonal mxm matrices with integer 
entries, so we can write in components 

t a(3 = Qa^aP («, P = 1 , • • • , TH, no SUH1 Over a) . (2.1) 

Introducing a U(l) gauge fields A® for each of the generators we define the covariant 
derivative 

m 

Df,w = d„w + i^2t a A*w, (2.2) 

a=l 

and the curvatures 

F% = d^A a v - d u Al. (2.3) 

Here the space-time indices v, .. run over {0, 1, 2}. In the following we shall exclusively 
work on (2+l)-dimensional Minkowski space with signature (—,+,+), whose points we 
denote by x M . Occasionally we also use latin indices i,j,k... 6 {1,2} to label the spatial 
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components of x M , or polar coordinates (p, 9) for (xi,x 2 ). Finally we should define a 
suitable "Higgs potential" for each generator t a of the gauge group U(l) m . The appropriate 
potential turns out to be the Fayet-Iliopoulos D-term (R a — wH a w) 2 ||, where R a is a 
parameter of dimension mass. There are a number of ways to think about this term, 
but we only point out as an aside that in terms of the symplectic geometry of C m , with 
symplectic form dw< A dw, H a — R a — wH a w is a Hamiltonian for the U(l) action on C m 
generated by t a . The interested reader is referred to the book || for more details on this 
point of view. 

Thus we can write down the general Lagrangian density which is the main subject of 
this paper. In natural units c = H = 1 it reads 

-i m -1 m 2 

C = --{D lt w)iD»w - ^( F P 2 - E Y (Ra " wHaW)2 - (2 ' 4) 

a=l a a=l 

It depends on m coupling constants e a of dimension (mass)^, and on the m parameters 
R a - The dependence of the theory on these parameters is one of our main interests. Note 
in particular that we can trivially eliminate any of the gauge fields by setting them 
and the corresponding coupling constants e a to zero. 

We are mostly interested in static fields, for which the energy functional has the form 
E = - d 2 x {D x w)^D lW + {D 2 w)^D 2 w + £ -(i^) 2 + ^ ^(R a - wh a wf . (2.5) 

a=l Ca a=l 

To ensure that the energy of a configuration is finite we also impose the boundary conditions 

lim (R a - wh a w) = (2.6) 

p — >oo 

lim D lW = 0. (2.7) 

p — >oo 

Below we will state more precisely how quickly these limits should be attained. For now 
we proceed, assuming that the decay of the fields at infinity is fast enough to justify the 
following manipulations. Thus using the algebraic identity 

1 f m 
E=- d 2 x\(D 1 ±iD 2 )w\ 2 + J2\- F i2±^-(Ra-wH a w)\ 2 

a_1 (2.8) 
If m 

/ d 2 xi(D lW y D 2 w - i(D 2 w) ] D lW + ^F 1 a 2 (i? a - wh a w) 

J a=l 



and integrating by parts we deduce 

1 f m 1 

E=- d 2 x\{D 1 ±iD 2 )w\ 2 + J2\-F? 2 ±^(R a 

ZJ a=l ea 1 

1 m 

4i> 



2 a- , 



where we have also used that 



(-D1-D2 - D 2 D 1 )w = 1 wh a w F± 2 



a=l 



Finally defining suitably normalised magnetic fluxes 



and 



we deduce the inequality 



More precisely we have 



o=l 



E > n\T\. 



E = ttT 

if and only if the Bogomol'nyi equations hold 



(D 1 +iD 2 )w = 

e 2 

F? 2 + -±(R a -wH a w) = 0, 



and 

E = -ttT 

if and only if the (anti)-Bogomol'nyi equations hold 



(Di -iD 2 )w = 

^I a 2-y(^a-^a^)=0. 



The quantity T deserves two further comments. First note that, by converting the 
surface integral into a line integral round the circle C at spatial infinity and using the 
boundary conditions ( |2.6|) and ( |2.7|) , T can be written as 

T = — (p wUw. (2.18) 
2m J c 

The second comment is a caveat. Since we are working on non-compact R 2 the fluxes 
$ a are not necessarily integers. Thus neither the fluxes nor T have, in general, a clear 
topological meaning. Nonetheless the magnetic fluxes are interesting quantities to consider 
because they are conserved if one rules out infinite energy configurations. This follows from 
Faraday's law of induction 

where Ef = is the electric field of the a-th gauge field. The integral on the right hand 
side is only non-zero if the electric field falls off for large p no faster than 1/p, which is 
precisely the condition for the electric field to have infinite energy. 



3. The Bogomol'nyi equations on R 



For the rest of this paper we focus on the Bogomol'nyi equations ( |2.15a|) and ( |2.15fc|) . 



We are interested in finite energy solutions of these equations, and we need to to specify 
more carefully the boundary condition which ensure that the energy and flux integrals 
written down in the previous section converge. For this purpose it is convenient to introduce 
the complex notation z = X\ + ZX2, d z = \{d\ — 182) and A a = \{A1 — iA%). Further we 
define gauge potentials and curvatures 



(3.1) 



Q>a ^ Q aa.A a 

a=l 
rn 

fl2 = QaaF 12 
a=l 

and we introduce the set of parameters r a through 

m 

Ra ^ ^ Qaa^ct (3.2) 



a=l 
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(note that the matrix Q with entries Q aa is invertible by virtue of the generators t a being 
independent) . The parameters r a are convenient for discussing the spontaneous symmetry 
breaking in this model. The potential 



m 9 

e z 



Y^^(R a -wh a w) 2 (3.3) 

vanishes if and only if wH a w = R a for all a or, equivalently if 

\w a \ 2 = r Q , a = 1, ra. (3.4) 

Thus, in order to have any chance of finding finite energy solutions we must require 

r a > 0, a = l,...,m. (3-5) 

Further we see that the U(l) gauge group which rotates the phase of wp for some given 
(3 G {1, m} is spontaneously broken if rp > but unbroken if rp = 0. 

In the new notation the Bogomol'nyi equations ( |2.15a| ) and ( |2.156| ) become 

(dz + ia a )w a = (3.6a) 

^ m 

fi2 + g E e lQ^Qap(rp - \wp\ 2 ) = 0. (3.66) 

a, 13=1 

These two first order equations imply one second order equation for w a as follows. At 
points where w a ^ the first of the above equations is equivalent to 

a a = idzhiw a , (3.7) 

which in turn implies 

/f 2 = iAln|«; a | 2 . (3.8) 
To extend this equation to the whole plane one uses the <9-Poincare lemma in a standard 



fashion, see e.g. |]10[| . Supposing that w a has zeros at {z ayS \s = 1, ...,n a } and using that 
in two dimensions 

Aln|z - z s \ 2 = 4n5(z - z s ) (3.9) 

one concludes 

A a 2 = -Aln|«; a | 2 - 2tvJ2^ ~ z a , s ). (3.10) 

s=l 
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Combining this with (|3.6fc| ) we arrive at the promised second order equation for w a : 

Aln|u; Q | 2 + e 2 aQa a Q a p(rf3- \w p \ 2 ) =47r^5(z-z Q)S ). (3.11) 

a,/3=l s=l 

This equation should be thought of as a generalised vortex equation, in the sense that 
in the case m = 1 and t\ = 1 it reduces to the equation 

N 

AlnH 2 + e 2 (r- \w\ 2 ) = 4tiJ2 6 ( z ~ z °) ( 3 - 12 ) 

s = l 

(we have omitted the label a = 1 and written iV = ni), which is the key to the study 
of Nielsen-Olesen vortices in the abelian Higgs model at critical coupling. As mentioned 
in the introduction many of the questions we will address in this paper are motivated by 
what is known about the abelian Higgs model; it may thus be useful briefly to summarise 
some salient results. 

The first question that one needs to settle is that of existence and uniqueness of 
solutions. In the case of ( |3.12|) this question is completely answered in ||1Q|| . It is shown 
there that no finite-energy solutions exist if r < 0, but that for r > 0, ( |3.12| ) has a unique 
finite energy solution for each integer N > and given zeros z s , s = 1, ...,N. It follows 
that for given N > (and thus fixed energy) the Bogomol'nyi equations in the abelian 
Higgs model have an iV complex parameter family of solutions, the iV complex parameters 
determining the positions of the zeros of the w-field up to permutations. This parameter 
space is called the moduli space of iV-vortices. The moduli space is in fact a differentiable 
manifold and has a natural Riemannian metric, inherited from the field theory kinetic 
energy functional. This metric is crucial in the so-called moduli space to vortex dynamics. 
In that scheme the interacting dynamics of iV vortices is approximated by geodesic motion 
on the moduli space of iV-vortices. In [|J it is shown that the Bogomol'nyi property of the 
field theory implies that the metric on the moduli space is Kahler. The Kahler property 
together with symmetries is then used to compute the metric explicitly in the case N = 2; 
by studying the geodesies of this metric much can be learnt about the interaction of two 
vortices. The Kahler property is also crucial in the study of statistical mechanics of vortices 
on Riemann surfaces in ref. 0. 

The example of vortices in the abelian Higgs model shows paradigmatically how the 
Bogomol'nyi property is the key to understanding difficult and interesting aspects of vortex 
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physics, such as vortex interactions and statistical mechanics of vortices. The family of 
theories defined in this and the previous section also have soliton solutions of Bogomol'nyi 
type whose physical properties can be studied with similar methods. In this paper we will 
not attempt to do this in full generality but instead go through some of the steps outlined 
above for a few specific model. However, to end this general section we show how one 
can deduce some general information about the spectrum of Bogomol'nyi solitons directly 
from ( FTTH ). 

In particular we can now state the conditions for a solution of the Bogomol'nyi equa- 
tions to have finite energy. The energy of solutions of the equations ( |3.6q| ) and ( |3.6fc| ) can 
be written conveniently in terms of the fluxes 

<p a = - 7 L [ d 2 xf? 2 . (3.13) 



2tt . 
It is 

m 

E = nT = n^r a tp a . (3.14) 



a=l 

2\ 



It thus follows from ( |3.6fc| ) that all fluxes and the energy are well-defined if (r Q — \w a \ ) is 



integrable over R 2 for all a. If r a > for all a the gauge symmetry is completely broken. 
Then all scalar and gauge fields approach their vacuum values exponentially 

\w a \ 2 »r a + C a e~ maP for large p, (3.15) 

with m a > and C a arbitrary constants, and there are no convergence problems. On 
the other hand if rp = for some (at least one) (3 G {1, ...,m} then the corresponding 
scalar fields wp and (because of the coupling) possibly some other scalar fields approach 
their vacuum values slower than exponentially. The precise asymptotic behaviour depends 
on the matrix (XllILi e aQaaQ a /3) which appears in ( |3.11|) and which couples the various 



components of w; we will give a more detailed discussion of some special cases later in 
this paper. Here we note that it is consistent with ( |3.11| ) for \w a \ to approach the vacuum 
according to a power law 



\w a \ » r a + C a p 2ria for large p. (3.16) 

In that case we impose 

??a > 1 (3-17) 
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to ensure finiteness of the fluxes and the energy. Finally it is possible that the power- 
law decay is modified by logarithmic terms 0. In particular we should also allow for the 
asymptotic form 

C 

\w a \ 2 « r a + " for large p, (3.18) 

p A In p 



which (for suitable e a and Q aoi ) is consistent with (|3.11|) and the finite energy requirement. 

We can obtain explicit formulae for the fluxes from (|3.10| ) as follows. Using Stokes's 
theorem we first obtain 

f 9 

/ d 2 x A In \w a \ 2 = 2ir lim p—ln\w a \ 2 . (3.19) 
J p^oo dp 

Then using ln(r a + e) ~ In r a + e/r a if r a ^ and e small, and integrating (|3.10|) we deduce 



r a > => (p a = n a e Z^ u (3.20) 

regardless of how |u> a | 2 approaches r a . If on the other hand r a = we know that |u> a | 2 
tends to zero according to (|3.16|) or (|3.18|) . Then the result of the integration can be 



summarised in 

r a = =>■ (p a = n a + rj a , n a G Z-° and r) a G R- 1 . (3-21) 

The real number rj a equals f] a if \w a \ 2 decays according to the power law ( |3.16| ) and is 1 
if \w a \ 2 approaches zero according to ( |3.18[ ). 

The above formulae show that the flux ip a counts the zeros of w a with multiplicity. If 
r Q > 0, all zeros are at finite p and have integer multiplicity. If r a = then w a has a zero 
at infinity whose multiplicity is > 1 but not necessarily integer. (One shows similarly that 
for square-integrable solution of eqs. ( |2.17q| ) and ( |2.176| ), (p a < if r a > and (p a < — 1 
if r a = 0.) A further useful condition can be deduced from ( |3.66| ) if rp = for some 
f3 G {1, m}. It then follows that 



a=l 

so that 



E(Q _1 W^>° (3-22) 



E(Q _1 )^^^°- (3-23) 

o=l a 



Here the equality holds only if wp vanishes everywhere. 



1 This was pointed out to me by Trevor Samols 



10 



4. Semilocal vortices, CP 1 lumps and toric varieties 



For specific calculations we will mostly concentrate on the GLSM with two scalar 
fields w\ and w<i in this paper. Our favourite choice of generators of U(l) subgroups of 
the maximal torus in this case is 

(we will mostly omit writing the identity matrix t\ in the following formulae). Then, to 
avoid crowded notation, we write and for the gauge fields A 1 and A 2 A and B for 
the corresponding complex fields A 1 and A 2 , and and for the corresponding field 
strengths. 

In this section we restrict attention to the case where one of the gauge fields is set 
to zero. Specifically consider the case e2 = and B^ = 0. The resulting model has 
been much studied in the recent literature and has vortex solutions known as semilocal 



vortices, see refs. |4[], flTj] and ||12|1 . One interesting aspect of this model is that it contains 



stable vortices although the vacuum manifold, defined as the submanifold of C 2 where 
the potential {R\ — w^w) vanishes, is a three-sphere (provided Ri > 0) and therefore has 
trivial first homotopy group. The reason why there are nonetheless stable vortex solution, 
first explained in [|J , is that in addition to condition ( |2.6| ) we have the condition (|2.7|) and 
this forces w to lie both on the vacuum manifold and on a gauge orbit of the gauged U{1) 
at spatial infinity. Thus, in terms of spatial polar coordinates (p, 6) w has to be of the 
form (wi, W2) = (cie lNl ®, C2e %Nl9 ) for large p with N\ some integer and c\ and C2 complex 
constants satisfying c 2 + c 2 , = R\. The point is that the gauge group acts without fixed 
points on the vacuum manifold, so that gauge orbits are necessarily loops and never just 
a point. The integer Ni is the degree of the map w\ p=00 from the circle at spatial infinity 
into one of these gauge orbits. Evaluating the formula ( |2.18| ) for the topological lower 
bound in this case we find T = R1N1, showing that vortex solutions of the Bogomol'nyi 
equations with N\ 7^ cannot decay into the vacuum. 



It has also been observed in the literature [0 that there is a close connection between 
semilocal vortices and topological solitons in the CP 1 (or 0(3)) sigma model. This con- 
nection is usually discussed in terms of energy scales, with the sigma model being thought 
of as a low-energy effective theory of the vortex model. However, from the present point 
of view it is more convenient to keep the energy under consideration fixed and vary the 
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parameter e\ (recall that this has dimensions (mass) 5). In fact it is easy to see that in 
limit e\ — > oo the semilocal vortex model reduces to the CP 1 model. To keep the energy 
finite we simultaneously impose the constraint 

R 1 -w^w = 0. (4.2) 

Then, since the kinetic term of the gauge field disappears form the Lagrangian in the 
limit e\ — > oo we can eliminate the gauge field altogether form the energy functional via 
the equation 

BC 1 

-((D lx w)^w-w^D fX w)=0, (4.3) 



dA^ 2i 



which implies 

R x Ap = iw^d^w. (4.4) 



It is a standard result, see e.g. |14]], that the resulting model is equivalent to the CP 1 sigma 
model. Geometrically, the condition Q4.2|) forces w to lie on a three-sphere S 3 , and the 



appearance of the covariant derivative with the gauge potential given by ( |4.4j ) means that 
the Lagrangian depends on w only up to an overall phase. Thus, defining the equivalence 
relation ~ via (wi, w^) ~ e lx (wi, w^), x G [0, 27r), the limit e\ — > oo leads to the non-linear 
sigma model with target space (S 3 / ~) = CP 1 = S 2 . For later use, we note that this can 
be made explicit by introducing the CP 1 - valued field 

u = — (4.5) 
or the 5" 2 - valued field cf> = (0i, fa, 03 ), defined via 

<Pi=w^nw, 1 = 1,2,3, (4.6) 
where t\,t% and Ts are the three Pauli matrices. The Lagrangian density can then be 



written in terms of u or 6, see ref. 14 



Note also that the topological character of the topological defect changes as we take 
the limit e\ — ► oo. For the semilocal vortex we have already interpreted the topological 
bound T as [K\ times) the degree of the w viewed as a mapping from the circle at spatial 
infinity to a gauge orbit on the vacuum manifold. However, in the limit e\ — ► oo, points on 
gauge orbits should be identified, so that now w maps the entire circle at spatial infinity 
into on point on the target space CP 1 . This allows us to regard « = W2/W1 as a map 
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from compactified space R 2 U {00} = S 2 to CP 1 . The flux $x = -1/2tt / d 2 xF 12: with 
F12 computed from ( |4.4|) , equals the degree of that map. Thus, in the general terminology 
of the introduction topological solitons in the CP 1 model are textures. However, in this 
particular context they are more usually called lumps. 

The moduli spaces of both semilocal vortices and CP 1 lumps have been studied in 
some detail. In [11] it is shown that the moduli space of semilocal vortices with magnetic 
flux number N is diffeomorphic to the space of polynomials of the form Pn{z) = z N + 
^2^ = q a s z s and Qn(z) = ^2^=o b s z s . The coordinatisation of that space in terms of the 

27V 

complex coefficients a s and b s , s = 0, N — 1 shows that it can be identified with C 
The translation from a point in the moduli space to an actual field configuration, however, 
is best done in terms of the zeros of Pn and Qn- They are also the zeros of the scalar 
fields wi and w^- 

In the CP 1 model, a lump of degree iV is a rational function on R 2 of degree N which 
tend to zero at infinity. Explicitly such a function can again be written in terms of the 
polynomials Pn and Qn 

«M = ffy (4-7) 

but now we have to require in addition that Pn and Qn have no common zeros, or 
equivalently that the resultant of Pjv and Qn is non- vanishing. Writing Rn for the set 
where the resultant does vanish we conclude that the moduli space of degree N lumps is 

27V 1 

C — Rn- The interpretation of the moduli is now quite different from the vortex case. 
Writing for example a single lump configuration as b/ (z — a), the complex number a is the 
lump's position in R 2 and the modulus and phase of the complex number b are its size 
and orientation respectively. Note that b has to be non-vanishing and that in the limit 
I fc» I — > the lump becomes infinitely spiky. These features generalise to lumps of degree N. 
The moduli encode information about internal as well as position degrees of freedom. The 
condition of the non-vanishing resultant removes precisely those points from the moduli 
space which correspond to infinitely spiky configurations. 

As mentioned earlier, the Riemannian metric which the moduli spaces inherit from 
the field theory kinetic energy is crucial in the moduli space approximation to soliton 
dynamics. In ref. |TH] this approximation is applied to lumps in the CP 1 model and it 
turns out that the moduli space metric has two problematic features in this case (see also 
16(1 for further discussion of this point). The metric is not finite because changes in the 
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coefficient q^-x require infinite kinetic energy. Furthermore it is not complete: there are 
geodesies which reach infinitely spiky configurations in finite time. By contrast the metric 



on the moduli space for semilocal vortices , studied in ref. |12[ in the case e\ = 1, is finite 
and complete. Thus the moduli space of semilocal vortices can be thought of in a very 
precise sense as a regularised version of the moduli space of CP 1 lumps, to which it tends 
in the limit e\ —>■ oo. 

Let us briefly consider the case where only the L r (l)-factor generated by t 2 is gauged. 
Thus we set e± = and = 0. The Bogomol'nyi equations ( p.l5a|) and (|2.155|) are this 

case: 

(dz + iB) Wl = 

(d- z - iB)w 2 = (4. 



G 12 + l(R2-\w 1 \ 2 + \w 2 \ 2 ) =0. (4.86) 

The energy of a solution of these equations is E = R 2 Q 2 . Thus, if R 2 < and E ^ we 
necessarily have $2 < 0. However, it then follows from a standard vanishing theorem (a 
line bundle of negative degree cannot have a non-zero holomorphic section) that w\ = 0. 
Similarly if i?2 > we deduce that for any solution with non- vanishing energy w 2 = In 
either case the vortex solutions of this model are just embedded Nielsen-Olesen vortices. 
In the limit e 2 — > 00 we again obtain a non-linear sigma model, but this time the target 
space is {{w u w 2 ) e C 2 \R 2 = \wi\ 2 — \w 2 \ 2 }/ ~, where ~ is the equivalence relation 
(wi, w 2 ) ~ (e l %i, e~ lx w 2 ). This space is known in mathematics as the weighted projective 
space CP( X _ 1 ) and is isomorphic to C. Bogomol'nyi solitons in the corresponding non- 
linear sigma model are holomorphic maps from S 2 to CP^ _ x y However, such maps are 
necessarily constant maps. Thus there are no non-trivial textures of Bogomol'nyi type in 
this model. 

The phenomena encountered in this section form part of a very general story, much 
discussed in mathematics and string theory, see refs. § and ||. Briefly, it goes as follows. 
Returning to the general notation of sect. 2, consider taking the limit 

e a — > 00, a e I (4.9) 

and simultaneously imposing 

(R a - wh a w) =0, ael (4.10) 
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for some subset of indices I C {1, ...,m} which we can without loss of generality take to 
be / = {1, <i}, d < m. Assume that the R a are such that the entire gauge symmetry is 
spontaneously broken. Also let us at first restrict our attention to the situation where the 
couplings and gauge fields labelled by the complementary indices are set to zero 

e a = 0, A a ^ = for a = d+l,...,m. (4.11) 

Defining the currents 

51 = ^ {(D„w)U a w - wH a D^w) , (4.12) 
the equations of motion for the gauge fields A* aGl are then 

j* = for ael (4.13) 

which one can solve explicitly for A®, a E I: 

d 

wh a t h w A\ = iwUad^w. (4.14) 

6=1 

The gauge fields A^, a E /, can thus be eliminated, but the Lagrangian (|2.4| ) is still 
invariant under transformations w(x) — > e %x ^'°to(i), a E I, Xa G [0, 2tc). Then 
defining the equivalence relation 

w~ e lxata w, a El, Xa G [0, 2tt) (4.15) 

the fields w may be thought of as taking values in the non-linear space 

Z = {wE C m \wH a w = R a ,aE I}/ ~ . (4.16) 



The combined operation of imposing ( ^.10|) and dividing by the action (|4.15| ) of 
the torus U(l) d is called the symplectic quotient of C m by U(l) d and is often written 
C m / /U(l) d . In fact what we are looking at here is a very special sort of symplectic quo- 
tient. With the original space being C m and the group action being a C/(l) d -action the 
resulting quotient is a so-called toric variety of complex dimension m — d. The complex 
projective space and the weighted complex projective space which we encountered above 
are special examples of toric varieties. Toric varieties are naturally Kahler manifolds, and 
their topology (which depends on the values of the R a ) is well-studied. In particular it is 
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known that under some restriction on the ( |2.1|) the quotient space Z is compact. This 
holds for example if for some a G I all the Q a a are positive. (For more general conditions see 
0.) Moreover, the second homology group H 2 (Z, R) is <i-dimensional for generic values 
of R a (note that as long as d < m, the real dimension of Z is > 2). 

In analogy to our discussion of the CP 1 sigma model static fields which obey the 
condition (|2.7|) can here be regarded as maps from R 2 U {00} = S 2 into Z. Such maps 
can be classified topologically by their multi-degree, the integral of the pull-back of the d 
generators of H 2 (Z, R). This generalises the observation made in the simple case of the 
CP 1 sigma model. For every gauge field which we eliminate by taking the corresponding 
coupling constant to infinity the topological meaning of the magnetic flux number changes. 
If before elimination this number counts the number of times the circle at spatial infinity 
is wrapped round a certain loop in the vacuum manifold then, after elimination, it counts 
the number of times compactified space get wrapped around a certain generator of the 
second homology of the (now non-linear) target space. 

Note that the argument in the last paragraph depends crucially on the fact that all 
gauge fields are either eliminated through the limit ( |4.9| ) or set to zero ( |4.11| ) . Only in this 
case does the boundary condition ( |2.7| ) allow us to identify points at infinity to one point. 
However, it is also interesting to consider the mixed situation, where some gauge fields are 
eliminated and others remain as dynamical fields. This leads to gauged non-linear sigma 
models, where the geometric interpretation of the magnetic fluxes f a is more subtle. We 
shall see this in a particular example in the next section. 



5. Topological solitons in U{1) x U{1) gauge theory 

We now have all the ingredients necessary for the study of the case where the entire 
maximal torus of U(2) is gauged. There are eight independent parameters in this model: 
the real numbers ex, ei, R\, R2 and the integer entries in the generators t a : 

It is instructive to write out the Lagrangian density in components for this model 



1,„ ,0 1,„ ,0 1 „o 1 

(5.2) 



'1 

- %R X - QnKI 2 - Q12M 2 ) 2 - ^(i?2 - Q21KI 2 - 

o o 
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2 „2 

2\2 



This Lagrangian density is of the type studied first by Witten in the seminal paper [17 



and later by many others (see the book |T8[ and references therein) as a model for bosonic 
superconducting vortices. In that context the discussion is usually conducted in terms of 
the gauge fields = QuA^ + Q2\B^ and = Q^A^ + Q22B^ which couple only to the 
scalar fields w\ and W2 respectively. Writing U(l) a and U{l)b for the gauge groups which 
rotate the phases of wi and W2 respectively the basic observation can be stated as follows. 
By adjusting parameters suitably one can arrange for one of the two gauge groups, say 
U(l) a , to be spontaneously broken while U(l)b remains unbroken. Thus U(l)b may be 
interpreted as the gauge group of electromagnetism (there is a caveat which we explain 
below). Then there is a range of parameters for which the model has £7(l) a -vortex solutions 
in whose core the Higgs field W2 of U(l)b is non-zero, thereby breaking the electromagnetic 
gauge group there. Witten showed that under these circumstances the core of the vortex 
becomes superconducting. 



The Lagrangian density ( |5.2|) is more general than the Lagrangians usually studied 
in the context of superconducting vortices because it allows for interactions between the 
gauge fields and (through the curvature terms) but it also has many special properties 
because it of Bogomol'nyi type. One general feature of Bogomol'nyi solitons in gauge 
theories with scalar fields is the cancellation between the static forces mediated by the 
gauge fields and those mediated by the scalar fields. Such a cancellation is of course only 
possible if the scalar and gauge fields are either both massless or both massive. Thus 
we should expect on general grounds that in our model, too, the scalar field W2 should 
acquire a long range component when U(l)b is unbroken. In fact it is easy to check this 
explicitly. By the general condition stated after eq. ( |3.5|) , U{l)b is unbroken if and only 
if T2 = 0, which is equivalent to Q21R1 = Qi\R2- The vacuum expectation value of 
\wi\ 2 is then Ri/Qn = R2/Q21 and, by collecting the terms quadratic in |w 2 |, the mass 
of W2 is found to be zero. (A similar calculation shows that if we fix the parameters so 
that w± vanishes in the vacuum then w± becomes massless.) From the point of view of 
superconducting vortex physics this result is disappointing: one requires an exponentially 
localised W2 condensate in order to interpret the unbroken gauge group U(l)b as the gauge 
group of electromagnetism. If the condensate only decays according to some power law 
the U(l)b gauge invariance is never properly restored outside the core of the vortex and 
the electromagnetic interpretation is not appropriate. 



By adding a suitable perturbation to the Lagrangian (|5.2|) one could ensure an ex- 
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ponentially localised condensate, of course at the expense of destroying the Bogomol'nyi 
property. However, properties such as interactive dynamics and thermodynamics of such 
'almost Bogomol'nyi' models can also be studied with relative ease by perturbation meth- 
ods ||19|1 . This may be interesting from the point of view of the phenomenology of su- 
perconducting vortices but we will not pursue it here. Instead we want to exhibit some 
detailed properties of the model ( |5.2| ) by picking a particular set of generators t a . 

Thus we return to the generators (|4.1|) , i.e. we set Qii = Q12 = 1 and Q21 = —Q22 = 1 
so that ri = (Ri + R 2 )/2, r 2 = (-Ri — -R2V2. We also continue with the index saving 
terminology and write a^ = Ai + Bi and hi = Ai — Bi, a = (a± — ia^)!^ and b = {pi — 7&2)/2 
as well as f\2 = F12+G12 and g\2 = F12 — G12. Then the Bogomol'nyi equations (3.6 )read 

(dz + ia)wi = (5.3a) 
(d z + ib)w 2 = (5.36) 
fi2 + \{e\ + el)(n - hi| 2 ) + ^(e? - e 2 2 )(r 2 - \w 2 \ 2 ) = (5.3c) 

912 + \{e\ - e\){r x - \ Wl \ 2 ) + ^(e\ + e 2 2 ){r 2 - \w 2 \ 2 ) = 0, (5.3d) 

and the generalised vortex equations (|3.11| ) are 

AlnKI 2 + (el + ej)(ri - K| 2 ) + (e? - e 2 2 ){r 2 - \w 2 \ 2 ) = 4tt ^ S(z - z hs ) (5.4a) 

s=l 

"2 

Aln |«; 2 | 2 + (e 2 - e 2 )(n - \ Wl \ 2 ) + (e 2 + e 2 )(r 2 - |«; 2 | 2 ) = 4tt ^ 8{z - z 2ft ). (5.46) 

t=i 

We are interested in the way the properties of solutions, such as energies and magnetic 
fluxes, depend on the coupling constants e± and e2 and on the values of the parameters 
ri and r 2 . The latter determine the symmetry pattern. If n > and r 2 > the gauge 
symmetry is completely broken. Then the fluxes <pi and <p 2 equal the integers ni and n 2 
respectively. Although we do not prove this rigorously here we expect that given these 
integers there is a (ni + 112) complex parameter family of solutions of (|5.4a| ) and ( |5.46| ). 
Each solution is fully characterised by the unordered set of parameters {^1,1, — , ^i, ni } and 
{z2,i, Z2,n 2 } which label the zeros of the scalar fields w\ and W2- Note that, contrary to 
the situation for a single vortex, the magnitudes of the magnetic fields are not necessarily 
maximal at these zeros. 

If ri > and r 2 = 0, U(l) a is broken but U(l)b remains unbroken (and vice-versa for 
ri = and r 2 > 0). Thus ifi equals again the integer n\ which counts the zeros of wi, but 
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f2 may have a non-integral part. The interpretation, explained in sect. 3, is that W2 has 
a zero at infinity whose order is at least one but may be non-integral. Thus only [f>2 — 1] 
zeros (where [992] denotes the integral part of ^2) can be placed arbitrarily in R 2 U {00}. 
For given ipi = m and tf>2 we therefore expect there to be a n\ + [992] — 1 complex parameter 
family of solutions of (|5.4a|) and Q5.46Q . 



Finally in the special case r\ = t2 = the full gauge symmetry survives, but we 
shall see that the model has no finite energy solutions in this case. Before we discuss the 
spectrum of solitons in more detail we will now show that our model contains, in the limit 
ei — ► 00, the gauged 0(3) sigma model introduced in ref. M and analysed further in ref. 



5.1. The gauged 0(3) sigma model revisited 

In taking the limit e\ — > 00 we can now follow the general recipe of sect. 4. Thus 



we impose simultaneously the constraint |u>i| 2 + 1 2 = R\ > 0. Then we solve for A 



according to (4.14) 



R\A^ = iw^d^w — w^t^wB^. (5-5) 

Since the field (iwi, W2) is now on ly defined up to an overall phase we can again discuss the 
model in terms of the ratio u = W2/1V1 Q4.5| ). In particular the energy functional for static 
fields takes the form 

rn 1 f j2 1 ADxu\ 2 + \D 2 U\ 2 1 2 2 f r l\ u \ 2 ~ r 2\ 2 r Ka \ 

2=2] + (ITRF + 4 Gh + e K^r) ' (B ' 6) 

where the covariant derivative on u is DiU — (di — 2iB{)u. This functional is a generalised 
form of the energy functional of the gauged 0(3) sigma model introduced in ref. and 
formulated there in terms of the field (p (|4.6| ) . To recover the formulae in ref. one should 
set r% = 1, ri = and e2 = 1, and identify —2Bi and — 2O12 with what is called A4 and 
F12 there. Then the energy functional ( |5.6|) is 4 times the energy functional studied in ref. 
and the magnetic flux 

1 r ,o 



$2 = -— / d z xG 12 (5.7) 
2n J 

is -g- x the magnetic flux defined in ref. , which is denoted $ there. 
Again we have the inequality 

E > 7r(riv?i +r 2 (^ 2 ), (5.8) 
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and from eq. (5J3) one finds the following formulae for the fluxes 



1 f 2 DiuD 2 u - D 2 uDiu 1 , 2 — i^i-i 
fi = 7^— d x — / d x (5.9) 

ZITl J (1 + 27T J 1 + |up 

and 

1 f 2 DiuD 2 u - D 2 uD x u 1 /" 2 2Gi2 1fV . 
f2 = 7T-. d x + — d x ; |2 . (5.10) 

27TZ 7 (1 + |Mp) 2 27T 7 1 + |W| 

Such integral formulae are useful for explicit computations, but to understand the geomet- 
rical meaning of the magnetic fluxes (fx and <^>2 in this model it is best to recall how the 
magnetic fluxes are related to the numbers of zeros of the scalar fields in the vortex model. 
Since T\ and r 2 should not both be zero let us for definiteness assume that n > 0. Then 
(pi is an integer (called ri\ above) which counts with multiplicity the number of zeros of 
w\ or equivalently the number of poles of u. The flux (p 2 counts the zeros of w 2 and hence 
of u. If r 2 > 0, i.e. if the only remaining U(l) gauge symmetry is spontaneously broken, 
then |it| 2 = r 2 /ri at spatial infinity and all zeros of u are at finite p. Then f> 2 is also an 
integer (called n 2 above) which counts these zeros with multiplicity. If r 2 = 0, i.e. in the 
case of unbroken gauge symmetry, the condition ip 2 e R- 1 tells us that u must have a zero 
at infinity of order at least 1. 

As an aside we note that integer (pi is called degree of u in ref. @. When ip 2 is not 
an integer this is somewhat misleading because the scalar field u then decays at spatial 
infinity according to some power law with a non-integral exponent. Such a field cannot 
be viewed as a continuous map from R 2 U {oo} to CP 1 and therefore does not have a 
well-defined degree. As explained here, ipi should be thought of more generally as the 
number of poles of u, counted with multiplicity. 

The Bogomol'nyi equations in this model may be derived either by determining the 
condition for the equality in (|5.8|) to hold, or by subtracting ( |5.3a| ) from (|5.3fc| ) and ( |5.3d ) 
from ( |5.3af ). In either case the result is 

(d s -2iB)u = (5.11a) 
2 ( ^iM 2 - r 2 



Gi2+e H i+N 2 J =0 - (5-116) 

These equations imply a second order elliptic equation for u which can again be derived in 
two ways. Either, as first shown in ref. ||20|1 , from (|5.11a|) and ( |5.11fc| ) via the 9-Poincare 
lemma or by subtracting eq. ( |5.4a ) from eq. ( 5.46Q . In either case the result is 

(I I 2 — \ n l n 2 

ri ' l 7 | |2 r2 ) = -4tt Hz - + ^ Yl 5 ^ z ~ ( 5 - 12 ) 
' \ U \ ' s =l t=X 
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This equation with r2 = is analysed carefully in ref. [BD1. It is shown there that in 
that situation and for given fi and <p2 it has a ipi + [^2] — 1 complex parameter family 
of solutions, which agrees with our general counting argument given in the discussion of 
( |5.4a| ) and ( |5.45|) above. 

Physically one may think of a solution of the Bogomol'nyi equations fl5.11a| ) and 
( |5.116|) as a texture carrying magnetic flux 2$2 = fi ~ V2 which counts the difference 
between the number (with multiplicity) of poles and the number of zeros of u. In the 
case of spontaneously broken gauge symmetry the flux is quantised and can be positive, 
negative or zero. In the unbroken case, however, it is clear from ( |5.7| ) and ( |5.11fc| ) with 
T2 = that the magnetic flux $2 is positive for non-trivial solutions. Thus we deduce that 
$2 takes values in a finite interval 

< 2$ 2 = fx ~ V2 < fx ~ 1- (5.13) 

This formula was derived independently by Samols (as referred to in ref. 0) and by Yang 
in ref. |2(J . In the next section we shall see that it is a consequence of the general inequality 

At the end of the introduction I mentioned that the wish better to understand the 
mathematical structure and physical significance of the gauged 0(3) sigma model was the 
starting point of this paper. In this section we have seen that mathematically this model is 
a limiting case of a certain gauged linear sigma model with Fayet-Iliopoulos potential terms. 
Since the relevant gauged linear sigma model has solutions which describe superconducting 
vortices of Bogomol'nyi type this observation also sheds light on the physics. When r2 = 
the gauged 0(3) sigma model has an unbroken U(l) gauge group which we may identify 
with the gauge group of electromagnetism. In the core of a soliton solution the gauge 
symmetry is broken and in this sense we may think of the solitons in the gauged 0(3) sigma 
model as superconducting textures. As in our discussion of superconducting vortices in the 
previous section the Bogomol'nyi property of the model leads to a power-law localisation 
of the scalar condensate, which is phenomenologically disastrous . However, this problem 
can again be solved by moving away from the Bogomol'nyi limit. In fact a model of the 



required type is studied in ||21|| . The soliton solutions, called Skyrme-Maxwell solitons 
there, are exponentially localised. Asymptotically there is an unbroken U(l) gauge group 
which is broken in the centre of a soliton. Such a topological soliton may thus properly 
be called a superconducting texture. These textures only involve one scalar field and one 
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gauge field; from a mathematical point of view they therefore appear to provide a more 
economical model for superconducting topological defects than superconducting vortices. 

5.2. Symmetry breaking patterns and the soliton spectrum 

The purpose of this subsection is to highlight some consequences of the simple energy 
formula 

-Ebps = ttIVi^i + r 2 y? 2 ) ( 5 - 14 ) 

for Bogomol'nyi solitons in the m = 2 GLSM. As we saw above, this formula is valid for 
all values of e± and e 2 , including the non-linear sigma model limit e± —>■ oo (although 
the interpretation of (p± and <£> 2 changes in this limit). Again we organise the discussion 
according to the symmetry breaking pattern. 

First consider the case of completely broken symmetry, i.e. r\ > and r 2 > 0. Then 
fluxes fi and <^> 2 take arbitrary non- negative integer values n\ and n 2 . The resulting 
spectrum depends in an interesting way on whether n/r 2 is irrational or rational. In the 
former case there is no degeneracy in the energy level for given (p i and <£> 2 other than 
that coming from the arbitrary positions of the n\ and n 2 zeros of w\ and w 2 discussed 
earlier. When ri/r 2 is rational, however, solutions with different flux quantum numbers 
may be degenerate in energy. More precisely, writing ri/r 2 = p/q, where p and q are two 
integers which are relatively prime, we deduce from formula ( |5.14| ) that the solution with 
(^li V ? 2) = (^1,^2) has the same energy as the solution with (pi,ip 2 ) = {n\ + g,n 2 — p) 
(provided these integers are still positive). In the case e± = e 2 , where the equations ( |5.4a| ) 
and ( |5.4fc| ) decouple this degeneracy is easily understood: if the masses of a single wi-vortex 
and a single w 2 - vortex have a rational ratio, then we expect degeneracies in the energy 
levels of superpositions of wi- and w 2 vortices. However, for e\ 7^ e 2 , and particularly in 
the limit e± —>■ 00, this degeneracy is more remarkable. In this limit we interpret ^1 as the 
degree of u and 2$ 2 as the magnetic flux carried by the texture; in terms of these ( |5.14j ) 
reads 

-Ebps = 71" {{n + r 2 )degree[w] - 2r 2 $ 2 ) . (5.15) 

Thus the degeneracy for rational ri/r 2 means physically that in the gauged 0(3) sigma 
model changing the degree of a configuration is energetically equivalent to changing, by a 
suitable number of units, the magnetic flux it carries. 



22 



If r± > and r<i = then £7(l) a is broken and U(l)i remains unbroken; in that case 
<Pi is still quantised as a non-negative integer, but f2 may now have a non-integer part. 
The energy (|5.14| ) depends only on <f\ and is degenerate with respect to changes in if 2- 
Note, however, that the range of <f2 is restricted by the finite-energy condition if 2 > 1 and 



the further constraint coming from (|3.23|) for non- vanishing solutions. In this case it reads 

p 2 - p 2 

1 2 / r -i r'\ 

</>2 < -n-, o^l- (5.16) 

Together these conditions force if 2 to lie in a finite interval. In particular they may force 
f 2 to vanish for certain values of e\ and e2- This holds in particular for e± = e2 where 
this result, too, is easily understood. The equation ( |5.4fc|) for W2 is then just the abelian 
Higgs vortex equation (|3.12j ). As mentioned in sect. 3 this equation has no solution other 
than the trivial solution when r = 0. In the non-linear sigma model limit e\ — > oo the 
formula ( |5.16j ) yields the promised re-derivation of ( |5.13| ). It now reads <f2 < (fx, which we 
combine with the earlier condition if 2 £ R- 1 to 

l<<P2<<Pi- (5-17) 

This is equivalent to the inequality (|5.13| ). Note in particular that this condition cannot 
be satisfied if f± = 1 and that there is therefore no solution with precisely one pole in the 
gauged 0(3) sigma model. This result was first derived for the spherically symmetric case 



in ref. and proved in general in ref. [20|. 



The reverse situation, with r\ = and r2 > can be discussed in an analogous fashion, 
and the case r\ = r 2 = does not lead to any interesting solution: now the condition ( |3.23| ) 
implies, for non-trivial solutions, ( |5.16| ) as well as 



e 2 - e 2 

Vi < 2 I 2V 2 (5.18) 
ef + ei 



which is incompatible with ( |5.16| ). Thus the trivial solutions wi = W2 = is the only finite 
energy solution of the equations ( |5.4q| ) and ( |5.4E| ) in this case. 

6. Abelian Chern-Simons models of Bogomol'nyi type 

Consider the following Lagrangian density for a theory where the dynamics of all the 
gauge fields are governed by Chern-Simons terms: 

1 m 1 

C CS = --(D^D^w - t—A^AI^ - Vesica, w, R a ), (6.1) 

l a=l lKa 

23 



with m dimensionless coupling constants K a . In this section we want to determine the 
potentials Vcs such that the model has static soliton solutions of Bogomol'nyi type. 

The Chern-Simons term is independent of the metric, so it does not contribute to the 
energy-momentum tensor. As a result the energy functional reads 

E CS = \J d 2 x(D w) j D w + {D lW ) ] D lW + {D 2 w) ] D 2 w + V CS - (6.2) 

However, the energy does depend on the gauge fields through the covariant derivatives. 
In particular the time-component of the gauge fields is determined as a function of the 
scalar fields and the spatial components of the field strength by the Chern-Simons version 
of Gauss's law: 

3$ = -^i a 2 , (6.3) 

where j° is defined as in ( f4.12| ). This constraint means that, even for static fields we cannot 
set Aq = 0. Instead we find, for time-independent w a , 

m ^ 

wUahw A b = —F? 2 . (6.4) 

6=1 Ka 

Thus the energy functional Ecs can be expressed entirely in terms of spatial components 
of the gauge fields and the spatial part of the field strengths. To do this it is convenient to 
introduce the notation F for the column vector with components F^/ K a , V for the column 
vector with components K a (Ra — w't a w) and W for the m x m matrix with components 
W a b = w't a tbW, a, b = 1, .., m. We claim that by setting 



Vcs = ^V*WV = - K a (R a - wH a w) wH a t b w K b (R h - wh b w) 

o o ~ — 

a, 6=1 

we obtain an energy functional of Bogomol'nyi type: 

E cs = \J d 2 x{D lW ) ] D lW + (D 2 w) j D 2 w + F*W _1 F + ^V*WV 

= - / d 2 x\(D 1 ± iD 2 )w\ 2 + (F ± i\VV) t W" 1 (F ± -WV) 
2, J 2 2 



(6.5) 



(6.6) 



i m r 

ToJ2 Ra d 2 xF? 2 . 

a=l J 



Here the superscript t means transposition and we have again used the boundary condition 
(|2.7|) to integrate by parts and exploited the identity ( |2.10|) . As in sect. 2 we therefore 
deduce the inequality 

Ecs > n\T\ (6.7) 
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with equality if and only if one of the Bogomol'nyi equations holds 



(Di ±iD 2 )w = (6.8a) 
± - wH a tbW Kb(Rb — wHi,w) = 0. (6.86) 



pa - m 
£12 

a 6=1 



In the case m = 1 the general prescription we have discussed leads to the abelian 
Chern-Simons-Higgs model with a sextic potential, which has been much discussed in 
the literature. One reason why this model has attracted attention is that it contains 
both topological and non-topological solitons p2]j. The generalised models described here 



include, for m = 2, the choice of generators (|4.1|) and Ki = 0, the semilocal Chern-Simons 



model discussed in ref. ||23| . and have in general a variety of both topological and non- 
topological soliton solutions whose properties deserve to be studied further. However, here 
we now turn to another point of principal interest. 

Like in the case of GLSM's with the gauge field governed by a Maxwell term it is 
interesting to consider limits in which some of the gauge fields are eliminated, thus leading 
to a (possibly gauged) non-linear sigma model. Specifically we can eliminate the Chern- 
Simons terms for the gauge fields A®, where a runs over the subset of indices / = {1, d}, 
d < m, by taking the limit 

K a — > oo, a E I (6-9) 

and simultaneously imposing 

(R a - wh a w) =0, a el. (6.10) 



Then the equations of motion for the gauge fields A^, a G /, are again ( 4.1 3| ) , allowing 



us to express these gauge fields in terms of the scalar fields and the complementary set of 
gauge fields: 

d m 

^wh a t b wAl = iwh a (d^ + i hA^w. (6.11) 

6=1 b=d+l 

However, there is a subtlety here: in the Chern-Simons case, unlike in the Maxwell 
case, simply taking the value of the potential in the limit ( |6.9|) and (|6.10| ) as the potential 
for the remaining degrees of freedom does not appear to lead to a model with interesting 
Bogomol'nyi solitons. The quickest way to see this is naively to apply the limit ( |6.9| ) 
and (|6.10| ) to the Bogomol'nyi equation (|6.86| ). The equations labelled by a G / then 
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lead to purely algebraic constraints on the scalar fields which eliminate many solutions. 
Here we therefore follow a different path and define a new potential as follows. Let I = 
{d + 1, m} be the complementary index set to / and let F and V be the column vectors 
with components Ff 2 /^5 and fta(Ra — wHaw) respectively, where a £ I. Further write 
W _1 for the matrix with elements (W -1 )-^, a, 6 G /. Then the kinetic term for the gauge 
field in the Chern-Simons energy functional in the limit (ft .Up is 



F t W" 1 F (6.12) 



and thus the appropriate potential is 



V—^-vUw- 1 ) 'v. (6.13) 



cs 8 

With this choice the energy functional 

E~ = - / d 2 x (DxvrfDxw + {D 2 w) ] D 2 w + F*W _1 F + -V* (w _1 ) V (6.14) 
satisfies 

E~ > tt\T\ (6.15) 
with equality if and only if one of the Bogomol'nyi equations holds 

(D 1 ±iD 2 )w = (6.16a) 

™a & - 

b=d+l 

In the first of these equations the gauge potentials Af, a G I are again determined by 
The potential (gig) is, in general, of a more complicated form than the sextic poly- 



nomial that is familiar in Chern-Simons theories of Bogomol'nyi type. Thus it may be 
useful to illustrate the general discussion by considering again the case m = 2, with entire 
maximal torus of U(2) gauged. We will again use the index saving notation of sect. 5. 

With t\ and t 2 as in ( f4.1|) , the matrix W is 

W=( W l W wH * W ) (6.17) 
\w't 2 w w l w J 
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so that in the limit K\ — > oo and with Ri — w^w = 0, the matrix W 1 is just a number: 

= p2 f\. , 2 - (6.18) 

Rf — (wU 2 w) 2 



Thus the Cher n- Simons potential is 

.2 



^cs = ^s-(Ri ~ (wh 2 wf)(R 2 - wH 2 wf (6.19) 

OXLl 



and the Bogomol'nyi equations are 



(D 1 ±iD 2 )w = (6.20a) 

G 12 ± - (wh 2 w) 2 )(R 2 - wH 2 w) = 0. (6.206) 

ZR\ 



Formulating this model in terms of the field u (|4.5|) is again instructive. Using also again 
the parameters r± and r 2 defined in sect. 5, the potential (|6.19| ) takes the form 



V CS = 2^(ri + r 2 y 1 v /' 1 ^ . (6.21) 

(1 + |w| J 

Choosing the upper sign in the Bogomol'nyi equations they now read 

(d 2 - 2zS)w = (6.22a) 
G 12 + 44(n + r 2 ) H }^~" 2) = 0. (6.226) 

In fact these equations and their soliton solutions are studied with slightly different notation 
in ref. [53] and ref. for the case of unbroken £7(1) gauge symmetry r 2 = 0. In ref. 



the case of broken gauge symmetry, r 2 > is also considered. As often in Chern-Simons 
theories there are both topological and non-topological solitons in this model. All soliton 
solutions carry magnetic flux, and in the case of broken gauge symmetry the topological 
solitons may derive their topological stability from the (quantised) magnetic flux, from 
the number of poles of u (the "degree" ) or from a combination of both. For more details 
we refer the reader to the papers ^4|, |25| and [^6| and also to ref. [^7|, where a related 



model is discussed. However, it should be clear from this brief case study that the family 
of gauged linear sigma models with Chern-Simons terms defined in this section have an 
even richer spectrum of solitons than their Maxwell cousins. 



27 



7. Discussion and Outlook 



In (2+1) dimensions, many interesting topological solitons of Bogomol'nyi type can be 
studied in a unified way in terms of the gauged linear sigma models discussed in this paper. 
Thinking of the GLSM's with Maxwell term and Fayet-Iliopoulos potential terms as a 
generalisation of the abelian Higgs model and keeping in mind the phenomena encountered 
here, it would be interesting to address the following points in full generality. 



Can one extend the analysis of ref. |I(J and ref. |2(J to establish rigorous general 
existence and uniqueness theorems for the generalised vortex equation (|3.11| )? If in the 
general model discussed in sect. 3 the gauge symmetry is completely broken (i.e. if all 
r a > 0) then all the fluxes (p a are non-negative integers and the examples studied above 
suggest that for given fluxes there is a Yl™=i complex parameter family of solutions. 
If the gauge symmetry is only partially broken, say r Q > for 1 < a < d and r a = for 
d < a < m, then I conjecture that there is generically a Q^a=i + X^Ld+itv 9 " — 1]1 
complex parameter family of solutions for given fluxes. However, in this case the fluxes also 
have to satisfy the constraints ( |3.23j ). Here we have only shown that these constraints are 
necessary conditions for finite-energy solutions. It would be interesting to know whether 
they (or possibly some stricter version of them) are also sufficient for the existence of finite 
energy solutions. 

Having counted the solutions it is then interesting to look at the structure - differen- 
tiable and metric - of the moduli spaces of solutions. The examples studied here suggest 
that the moduli spaces of vortices on R are typically diffeomorphic to C for some D 
(determined by the above counting arguments) and have smooth metrics. However, in the 
limit e a — > oo for some a G {1, m}, where vortices turn into (possibly gauged) textures a 
number of interesting things happen. The interpretation of the moduli changes. Whereas 
vortex moduli simply characterise the (unordered) zeros of the scalar fields, the moduli 
for textures include internal as well as position parameters. Taking certain limits of the 
internal parameters corresponds to making the texture infinitely spiky, and these limiting 
values are therefore forbidden as moduli for textures. Thus the moduli spaces for textures 
are typically obtained from those of vortices by removing certain lower dimensional alge- 
braic submanifolds. Metrically, the moduli spaces of textures are also worse behaved than 
those of vortices. The example of the CP 1 lumps suggest that non-complete and non- finite 
metrics arise. In studying moduli spaces of textures it is therefore useful to bear in mind 
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that they may be thought of as singular limits of smooth vortex moduli spaces. In that 
way, singularities may be understood and, if desired, circumvented. 

Some of the questions raised so far are completely answered in the case where GLSM's 
are defined on compact Riemann surfaces. There the counting of solutions can be done 
using standard theorems in algebraic geometry, and much can be deduced about the topo- 
logical and (complex) differentiable structure of the moduli spaces from the key observation 
that they (like the target spaces of non-linear sigma models encountered in sect. 4, and for 
similar reasons) are toric varieties; see ref. for a pedagogical explanation of this state- 
ment and for further references. In particular the intersection ring of the moduli spaces 
is described in the literature, see again ||. This is of interest from the point of view of 
vortex physics because knowing the intersection ring of the moduli spaces is the crucial 
ingredient in the study of the statistical mechanics of vortices in ref. . There the statis- 
tical mechanics of Nielsen-Olesen vortices is studied on (amongst other Riemann surfaces) 
the torus, which is equivalent to studying vortices on the plane with periodic boundary 
conditions. Thus it appears that much of the information needed similarly to study the 
statistical mechanics of the more general types of vortices described here can be found in 
the mathematical and string theory literature. However, one important ingredient in the 
analysis of ref. ||, a certain normalisation factor, is lacking. Computing it requires some 
knowledge of the Riemannian structure that the moduli space inherits from the kinetic 
energy of the field theory. Since that metric is only natural form a (2+1) -dimensional 
point of view it has not been considered in the string theory context. 

It is also interesting to ask some of the above questions for the Chern-Simons version 
of the GLSM's introduced in this paper. Already the first step, the general analysis of the 
soliton spectrum, is more complicated in this case. There are now both topological and 
non-topological solitons and since the scalar potential no longer has a universal polynomial 
structure the finite energy condition needs to be analysed for each model separately. It is 
also natural to contemplate further generalisations, such as mixed models, where some of 
the U(l) gauge fields are governed by Maxwell terms and others by Chern-Simons terms. 
It should even be possible to construct general abelian GLSM's of Bogomol'nyi type where 
both a Maxwell and a Chern-Simons term is included for each gauge field. As shown 
first for the abelian Higgs model in ref. [|28| , later extended in ref. [29] and recently 



demonstrated in the context of the gauged 0(3) sigma model in refs. [2JJ] and ||26|| , this 
requires the introduction of additional neutral scalar fields. In such models the gauge 
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fields can get their masses from the Chern- Simons term ("topological mass") or through 



the Higgs mechanism. The analysis in ref. |H§ suggests that generalising this construction 
for abelian GLSM's will lead to models which can accommodate a truly bewildering array 
of massless, massive, topological and non-topological particles. Thus it seems that the 
framework of abelian GLSM's provides an extremely versatile kit for the construction of 
Bogomol'nyi type field theories in (2+1) dimensions. The Maxwell version mainly studied 
here also shows that this kit comes equipped with powerful tools for exploring the physics 
of these field theories. 
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